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ANALYSIS OF THE BACKSCATTER SPECTRUM 
IN AN IONOSPHERIC MODIFICATION EXPERIMENT 

by 

Hongjin Kim 

Institute for Plasma Research 
Stanford University 
Stanford, California 


ABSTRACT 

The purpose of this study is to compare predictions of the back- 
scatter spectrum, including effects of ionospheric inhomogeneity, with 
experimental observations of incoherent backscatter from an artificially 
heated region. Our calculations show that the strongest backscatter echo 
received is not, in fact, from the reflection level, but from a region 
some distance below (about 0.5 km for an experiment carried out at 
Arecibo), where the pump wave from a HF transmitter 100 kw) is below 
the threshold for parametric amplification. By taking the standing wave 
pattern of the pump into account properly, the present theory explains 
the asymmetry of the up-shifted and down-shifted plasma lines in the 
backscatter spectrum, and the several peaks typically observed in the 
region of the spectrum near the HF transmitter frequency. 
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I. INTRODUCTION 

From experimental results obtained at-Arecibo (Carlson, Gordon and 
Showen, 1972; Kantor, 1972), it is believed that the enhanced heating 
and anomalous absorption occurring in the ionosphere are mainly due to 
parametric decay instabilities; the HF ordinary (pump) wave propagating 
in the F-layer from a high power transmitter (<*-. 100 kW) decays into 
electron plasma and ion acoustic waves to enhance the plasma fluctuations 
(DuBois and Goldman, 1965; Nishikawa, 1968). Several analyses have been 
carried out of the saturation spectra in homogeneous (DuBois and Goldman, 
1972a and b; Valeo, Oberman and Perkins, 1972; Kruer and Valeo, 1972; 

Harker, 1972; Kuo and Fejer, 1972; Fejer and Kuo, 1972; Perkins and 
Valeo, 1973) or inhomogeneous (Arnush, Fried and Kennel, 1972a and b) 
plasmas, with the assumption that the pump field is well above the 
threshold for parametric amplification. This results in saturation of 
the electron plasma waves due to their nonlinear Landau damping by 
positive ions, or due to their absorption by electrons whose orbits 
are nonlinearly perturbed by the waves (Bezzerides and Weinstock, 1972). 

However, some essential features of the backscatter spectrum still 
remain to be explained. For example, Carlson et al . (1972) have obtained 
a pair of plasma lines in the spectrum enhanced by the pump wave 
(fg = 5.62 MHz). These plasma lines are shifted up and down from the 
incoherent scatter diagnostic beam (f. = 430 MHz) by f -f , where 
f (^ 4 kHz) is the frequency of the ion acoustic waves generated, and 
differ from each other in the intensities of the peaks. As shown in 
Fig. 1, Kantor (1972) has also observed several maxima and minima near 
the peaks of the spectrum. So far, these features of the spectrum have 
not been explained. 

The purpose of the present analysis is to take the inhomogeneity of 
the ionosphere into account and to compare our predictions with observations 
of incoherent scatter. In contrast to most of the previous analyses, 
our theory assumes that the pump is below the threshold. Section II is 
concerned with the pump wave. Its intensity is derived as a function 
of height. In Section III, the intensity of the parametrically enhanced 
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FIG. 1. Frequency spectrum obtained by incoherent 
backscatter at Arecibo (Kantor 1972) . 


electron plasma wave is obtained from the kinetic wave equation. 
Section IV is devoted to the incoherent scatter experiment performed 
at Arecibo, and frequency spectra are obtained for several cases. Our 
assumptions are discussed and predictions are compared to experimental 
data in Section V. 
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II . INTENSITY OF THE PUMP WAVE 


As the pump wave from the HF transmitter propagates vertically 
upward in the ionosphere, its group velocity decreases very sharply 
near, and reaches zero at, the reflection level, as shown in Fig. 2. 

Since more slowly moving waves spend a longer time in a location, and 
are thus able to exchange energies more effectively, we shall assume 
that energy from the pump wave is deposited in the ionosphere only at 
or very near its reflection level and we neglect the depletion of its 
intensity due to the linear or nonlinear loss in the region of interest. 
The present investigation is primarily concerned with heating at fre- 
quencies below the penetration frequency, so that the electron density 
profile can be approximated by a linear profile, as shown in Fig. 2(a). 
For simplicity, we neglect the geomagnetic field. After the pump 
field is found, however, it is regarded as being aligned in the direction 
of the geomagnetic field. 

A vertically propagating ordinary wave, of angular frequency , 

is described by the wave equations 



€ ( z ) E = 0 





for z ^ 0 (in the ionosphere) 

( 1 ) 

for z < 0 (in the free space), 


where z is the height measured from the base of the ionosphere (z = 0), 
and is normalized to the scale height 11 q ; c is the speed of light; 
the wave magnetic field is related to the electric field by 


B 


laj o h o 


dE 

dz ’ 


( 2 ) 


and the plasma equivalent permittivity, e(z) , is given by 


:(z) = 1 - 


o> 2 

_R 

CD 2 

0 


= 1 - z 


(3) 
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FIG. 2. The electron density and the electron plasma frequency 
are shown in (a) and (b), respectively, as functions of z. 

(c) shows the synchronism parallogram for decay instability 
at different heights in the ionosphere. 







The solution to Eq. (l) yields 


E = a Ai (-£) + b Bi(-£) , for z ^ 0 , 

( 4 ) 

E = E 0 jexp(- i£^ /2 z)+R exp(i^ /2 z)J , for z < 0 , 


where Ai and Bi are the Airy integral functions (Abramowitz and 
Stegun, 1964); R is the reflection coefficient, and the variables Q 
and Cq a 1 "® defined by 


C = C 0 (i-z) , 


3 o 



(5.) 


The magnitude of the electric field intensity of the incident pump wave 
at z = 0 , Eq , is given by u 


E 


0 



(V/m) , 


( 6 ) 


where the average heater power is in watts ; h is the height of the 
base of the ionosphere from the ground, and the gain, G , of the antenna 
is given by 

f 2 0 , (7) 


where A is the aperture. In the absence of anomalous absorption, it 
is customary to put b = 0 in (4) because the electric field has to 
vanish as z -* ® (Ginzberg, 1964; Budden, 1961 ). In the present 
situation, where most of the pump wave energy is deposited at or near 
the reflection level, however, we assume that the electric field 
vanishes in z s 1 due to absorption, and phenomenogically keep the Bi 
term in the region below the reflection level (0 <■ z < 1). 

Application of the boundary conditions of continuity of E and B 
at z = 0 , and use of (2), reduce (4) to the form 
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2 


( 8 ) 


E = S E. 


O' 


where the swelling factor, S , is given by 


S = rtC 


1/2 


+ { exp 


| ex p[-i(f C 5 0 /2 - |)] +ReXP [ i (f " t).} Ai ( " C) 

[-(? ‘r • s)] 


+ R exp 


i[- r 5/2 + £ 

1 3 ^0 4 


Bi (-0 


(9) 


and use has been made of (Abramowitz and Stegun^ 1964) 


Ai(x)Bi / (x) - Ai'(x) Bi(x) = — 


( 10 ) 


Noting that £ » 1 , use has also been made of the following asymptotic 

expressions for the Airy integral functions valid for | x| » 1 


A.(-x) 


1/2 


1 -1/4 . / 2 3/2 it \ ' , \ 1 

x sm^— ~xr + , A^-x)**--^ 


1/4 

x cos 




, s l -1/4 

B (*x) « -jy 2 x COS . 

It ' 


(f 113/2 + 1 )’ B /-*> 


1/2 


1/4 . 

x sxn 


/ 2 3/2 it\ 

in (j X + 4)' 


( 11 ) 


The familiar results of S and R in the absence of the absorption 
are recovered from (9) by letting the Bi term equal zero (Ginzberg, 
1964; Budden, 1961), as 


s = *k)J 2 Ai 2 c-e> - 


R = exp[-i(| C 3 J 2 - I) 


( 12 ) 


At the reflection level (£ = 0) , • the swelling factor becomes 


O I 2.1/2 

s = 1 “ c o c o 


1+ R 


+ 2 1 r| cos I'— Cq / 2 + z + tj 


- 1/2 5 _2/5 

c Q = Ai(0) = 3 Bi(0) = 


(13) 


7 



where ^ is the phase of the reflection coefficient, R . The condition 
that the WKB solution be a good approximation to (1) is j^j» 1 . In 
this region, (9) becomes 

- - (^rj. ♦ |r | 2 + 2|n|coaj|(c5 /2 - e y2 )+*' 

Figure 2 shows how the swelling factor varies with height for small 
| R | , in the region where the WKB approximation is valid. 

The reflection coefficient, R can be measured from the echo of 
the pump wave. Hence, the field intensities of the pump wave at the 
reflection level, and in the region where the WKB approximation is valid, 
are determined by (8) with the swelling factors given by (l3) and (14-), 
respectively . 
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G. 5* Swelling factor for |R| = 1/2 and if = 0 as function of z. The frequency 
scales are obtained from ( 25 ). 






III. KINETIC WAVE EQUATION FOR THE ENHANCED PLASMA WAVE 


Figure 2(c) shows a set of electron plasma wave dispersion waves 

at different heights (z = 0 to l) and a locus satisfying the synchronism 

conditions for parametric interaction among the pump wave (u)_,k ), 

u -~u 


the electron plasma wave (d),k.) and the ion acoustic wave (d) ,k ). As 
shown in the figure, the electron plasma wave may be excited at any 
height in the ionosphere. We assume that the pump field is below 
threshold so that higher order processes such as the nonlinear ion 
Landau damping or electron orbit perturbations are negligible. The 
kinetic wave equation describing the excited electron plasma wave may 
then be written in a WKB approximation form as (DuBois and Goldman, 1972a; 
Valeo, Oberman and Perkins, 1972) 

2 

y l ukt 

( 15 ) 


■la) (It + vliK " (a?) 


“oY 2 e o^ 7, 


1+d 


0) 


& + 

k k D co 



1+d 


where v is the group velocity of the electron plasma waves; a is the 

—g 

square root of the ion electron mass ratio; k^ is the Debye wave number; 

(i is the direction cosine of k with respect to the pump field E^ which is 

aligned with the geomagnetic field; y and y are the linear damping 

' e ' a 

rates of the electron plasma and ion acoustic waves, respectively, and 



We define n^ as the electron density at the reflection height; © as 
the electron temperature in energy units, and 


X k = 


lim 
V,T - 


2m VT 


( 17 ) 


as the spectral density of the electron plasma waves. The dispersion 
relations for the plasma and ion acoustic waves are given approximately 
by 
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(18) 



where v is the electron thermal velocity. 

e _ 

The first and second (bracketed) terms on the RHS of (15) describe 
the parametric growth due to the pump wave and the linear damping, 
respectively. The third term describes the thermal fluctuations, which 
are small compared to the Cerenkov radiation term, the last one. Hence, 
we neglect the third term (Valeo, Oberman and Perkins, 1972). 

We now assume that the electron plasma wave has reached a constant 
amplitude at which the parametric growth is balanced by linear damping 
and energy loss due to propagation out of the interaction region. By 
using the expression for the group velocity 


v 

g 



1 



> 


Eq. (15) may be rewritten as 


Q 


[(%)' 


1/2 


aT = + Ps - 


■',lv 1+d 2 


, Q 


, 1/2 / ... 
2 » v e m 0 

2 7 HCU 


(19) 


( 20 ) 


where the direction cosine of k with respect to the z-axis, p', is 
related to p by 


p' = p cosBq - sin0 sin0Q cos cp , (21) 

where 0^ is the angle between the geomagnetic field and the vertical 
direction . 
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(1+d 2 )' 1 


in (20), the generated electron plasma wave has 


CD = ocu-cd 
0 a 


at any height . For a given o>, numerical solu- 


Since P 
its peak near 

tion for £^(z) is straightforward . The required boundary condition is 
e^(0) = 0 because.no electron plasma wave can be generated at the base of 
the ionosphere. 
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IV. THE INCOHERENT SCATTER SPECTRUM 


An incident diagnostic wave propagating vertically upward in the 
ionosphere is scattered by the enhanced electron plasma waves of ( 20 ) . 
The scattered wave detected by the receiver should satisfy the synchro- 
nism conditions (Bekefi, 19 66 ) 


co = 00. i (u 1 k 
SI ’ s 



up-shifted 

down-shifted 


( 22 ) 


For incident and scattered waves with frequencies very high compared to 

the penetration frequency of the ionosphere, to. is related to k. 

1 , s 1 , s 

by the simple dispersion relation 


co. = k. c . 

1, S 1,S 


( 23 ) 


The loci in the co, k-plane satisfying the synchronism conditions are ob- 
tained by eliminating k. and co from (22) and ( 23 ) as 

X y S S 

k s ± — | up-shifted J ( 2 4) 

(down-shifted) 


These loci are shown by heavy lines in Figure 4. Since co « 2co^, the 
last term on the RHS of (24) is usually dropped and the loci are considered 
as vertical lines at k = ±2k^ . However, it turns out to be an important 
term which causes the difference in magnitudes and shapes of up-shifted 
and down-shifted spectra. 

Figure 4 also shows two dispersion curves of electron plasma waves 
with angular frequencies co = co a -co for k given by (24) . Since, as 
shown in ( 20) , the electron plasma waves are enhanced most significantly 
near co = co^-co^ these two are among the most important waves. The 
figure clearly shows that they belong to different heights. The depen- 
dence of the frequency of electron plasma waves on height is found from 
( 18 ) and (24) as 


»(®) = 


/a> \ 2 12 “i'VV / V e\ 2 

(up-shifted ) 

(25) 

K) ' a) 2 U J 

( down-shifted ) 

> 
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FIG. 4. The synchronism parallelogram for up- and down-shifted 
plasma lines. 
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where, since a>. » cn. » cd , the value of k in (24) has been taken 
at ao = cDq. The heights far above two electron plasma waves are then 
found from (25) with oo = q u-oo : For typical values in the Arecibo 

experiment of f = 5*62 MHz, f^ = 430 MHz, f^ = 4 kHz, = 1200°K 
and H = 3 X 10^ m, they are 

(4.75) o (up-shifted ) 

h « H-J [ X 10 m , J . ( 26 ) 

(4.64' (down-shifted) 

The height where the maximum up-shifted plasma line comes from is about 
10 m below that for the maximum down-shifted plasma line. The magnitudes 
of the swelling factor at these heights are shown in Figure 3- Conse- 
quently, different intensities of the electron plasma waves at these 
heights are expected from (20), which in turn result in different magni- 
tudes of up-shifted and down-shifted spectra. The power scattered by 
the enhanced plasma waves is given by (Harker, 1972) as 


dP 

s 




PA 
t r 


(z+h/h Q ) 2 k 


dz , 


(27) 


where A^ is the receiver aperture area; P^ is the transmitted radar 
power, and the Thompson backscattering cross-section, a^, may be taken 
as the square of the classical electron radius for the present purpose, 


0 r = r o = 7-94 x 10 m 2 . ( 28 ) 

From (25); the variable z may be changed to go via 


CD 

dz = 2 —x dcu . 

aT 

0 

Equation (27) then yields 


dP 

s 

dcu 


2n.Q m v PA (2<b.± o>)o) 
0 T e t r v 1 ' 


°®Vo 


1/2 


< ? - h/ V J 




( 29 ) 


( 30 ) 
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where is n at the reflection height. The power spectrum is plotted 

in Figure 5 , with the phase of the reflection coefficient as a parameter, 
and with typical values for the Arecibo experiment: 

T = 1200°K , H = 30 km , h = 150 km , 

f Q = 5-62 MHz , 9 0 = 40° , 

f = Y /2it = k kHz , P. = 100 kW , 

3 a 0 (3 i) 

f = Y /2k = 65 O Hz , P_ = 2.5 MW , 

6 G 1 

f . = 430 MHz , R = 1/2 , 

Antenna diameter = 300 m , Antenna efficiency = 40$ . 
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FIG. 5- Frequency spectra (a) when \|r = 0 and n/2. 
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FIG. 5 (contd). (b) when ijr = it and 5^/2. 




V. DISCUSSION 


Our theory has been based on the subthreshold condition: in contrast 
to previous analyses, we have neglected nonlinear Landau damping of 
electron plasma waves by positive ions or effects of nonlinearly perturbed 
electron orbits as higher order processes. Since most of the energy 
deposition from the pump wave occurs at the reflection level (z = 1), the 
intensity of the enhanced electron plasma wave may be well above the 
instability threshold. However, as shown in (26), the strongest echo 
received does not come from the reflection level, but from a region about 
0.5 km below for typical experimental data obtained at Arecibo. In this 
region, the pump power from a 100 kW radar is about 0.8 times the 
threshold at most. In passing, we also note that the 'WKB approximation 
we have employed is valid in this region. 

Inclusion of the geomagnetic field may add to the RHS of (18) an 
2 2 

additional term m sin 0 , where u) is the electron cyclotron 

c c 

frequency due to the geomagnetic field. Though this is larger than 

the thermal correction term in the dispersion relation, assuming that 

0) sin 0 « 10 , and that the spatial dependence be much weaker than 

c 0 

that of the wave vector k , gives (20) and (25) as good approximations. 
Thus, the geomagnetic effect does not invalidate the present results. 

For a particular phase of the reflection coefficient, each spectrum 
in Figure 5 has peaks near f = ±( f -f ). These are due to electron 
plasma waves parametrically enhanced by the pump wave. The intensities 
of the up-shifted and down-shifted peaks differ from each other. This 
is due to the difference in heights where the scattered waves of f = 
±(f n -f ) are generated as calculated in (26), i.e., as shown in Figure 3; 
there is a difference in the swelling factors at these heights. While 
the down-shifted peaks are smaller than the up-shifted in Figure 5 for 
the phases of the reflection coefficient, = 0 and 3^/2, the spectra 
for = n/2 and jt, where the down-shifted peaks are larger than the up- 
shifted, correspond to the experimental observations. 

Consistent with the spectra obtained at Arecibo, there are several 
maxima and minima near the peak of each spectrum. These are due to the 
swelling factor. The frequency intervals between adjacent maxima ( 15 
~ IT kHz) are in good agreement with the experimental results ( 8 ~ 20 kHz) 
(Kantor, 1972 ) . 
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